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Abstract 


A Bose gas in an external potential is studied by means of the local density 
approximation. Analytical results are derived for the thermodynamic proper- 
ties of an ideal Bose gas in a generic power-law trapping potential, and their 
dependence on the mutual interaction of atoms in the case of a non-ideal Bose 


gas. 
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One of the most striking consequences of quantum statistics is the condensation of an 
ideal Bose-Einstein gas where the zero momentum state can become macroscopically oc- 
cupied at a sufficiently low temperature [fl]. Some of the features of the transition to a 
superfluid phase exhibited by the Bose fluid tHe is often interpreted essentially as a result of 
this Bose-Einstein condensation (BEC), the strongest degeneracy effect of a Boson system 
B]. For many years, it was considered hopeless to experimentally observe BEC in an atomic 
gas with weak interactions. With the development of techniques to trap and cool atoms, 
BEC was recently observed directly in dilute atomic vapors [BH]. The new experimental 
achievements have stimulated great interest in the theoretical study of inhomogeneous Bose 
gases. 

The thermodynamic properties of trapped atomic Bose gases undergoing BEC can be 
altered by the spatially varying trapping potential. The interaction between atoms may have 
a significant effect on the thermodynamic properties. There have been several investigations 
analyzing the dependence of the critical temperature on the trapping potential and weak 
interaction in the Bose gas [B]. The thermodynamic properties of Bose gases in an external 
potential have also been discussed in Ref. |f], but no analytic results have been given. Here 
we shall derive some analytical expressions for the thermodynamic properties of a trapped 
non-interacting and interacting Bose gas under the local density approximation. 

Let us first consider the case of an ideal Bose gas trapped in an external potential. When 
the energy level spacing due to the trapping potential is much smaller than the thermal 
energy kT = 67}, the local density approximation is adequate DHIT]. Space may then be 
divided into small cells, and in each cell we may consider the trapping potential V(r) to be 
constant. With a trivial modification to indicate r dependence explicitly we may directly 
extend the formula of the grand potential, which can be found in any textbook on statistical 
mechanics, e.g. Ref. [I], to write the local grand potential above the critical temperature T, 


as 


Opz(r) = — 5 95/2(C) (1) 
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where u is the chemical potential, the thermal wave length À = (27h? /mkT ) , and 


¢ =exp|-8(V(r)- u), a(x) =X j2. (2) 
j=l 
From this grand potential the local density of the Bose gas in the external potential is 
dQ px (r) 1 
(tn) = — (P0) Eol (3) 
This is the same as [[]]]. For the generic power-law potential discussed in Ref. 
x |P i z jl 
Veale te Paf (4) 


where Lı, Lə and L3 are the linear sizes of the volume, and €1, €2, €3, p, | and q the 
parameters, integrating the local grand potential ([]) over the whole volume, we find the 


total grand potential of the system to be 


Qe (Tp) = KTA” | gsja(G)ar 


W 
= ~ 3 gat t/2 gn+2(2), (5) 


where z = ef”, n = 1/p + 1/1 + 1/q + 1/2, and 


W = ae (p, l, q) = (pla)* P(1/p) P(1/t) P(L/a), (6) 
E1 Eq) E3 


In the derivation we have used the definition of the function g (x) and the formula for the 


gamma function 


T(z) = [ Ce dt. (7) 


From the total grand thermodynamic potential the average atom number can be obtained 


as 


da ðu 


W o 
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W 
= Sagar (2), (8) 
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This is in agreement with the result in Ref. [Q] obtained directly from the local density p(r) 
instead of the grand partition function. 


The total entropy can be obtained as 


o=- (GF) 


E EE [(n + 2) gn+2(z) — BHon+1(2)] - (9) 


Although, for obvious reasons, the volume and pressure are not useful thermodynamic vari- 
ables for a trapped gas |], we may still consider an analogue to the capacity for the ho- 
mogeneous system at constant volume with a fixed number of atoms, and define this heat 


capacity as 


OS 
C(T)=T (FF) (10) 


which includes the work done against the potential as the energy of the gas is changed. 

So far we have discussed the case for a temperature above T.. When the temperature T 
is below the BEC critical temperature T}, the chemical potential u = 0, and only the normal 
component of the Bose gas has a contribution to the heat capacity. From Eqs. (H) and (M0) 


we can then obtain the heat capacity (T < Te) to be: 


(7 +1)(n+2)Wk 


C-(T) = agar 9n42(1) 
= DUA ged (EY (N)k. (11) 


As a limit, at the critical temperature T,, the above expression reduces to 


(n + 1)(9 + 2) Gn+2(1) 


Oat 9n+1(1) 


Cc 


(N)k. (12) 


When the temperature T is above the critical temperature T}, at a fixed N the chemical 
potential u is a function of temperature T. In order to obtain the heat capacity, we have to 


calculate the derivative of u with respect to temperature T while holding N fixed. Using 


O(N) 


from Eq. (§) we find 


(n+ 1)gn+1 (z) + 


1 (Op 
a a nle) =0, (14) 


which yields 


(2%) _ k(t 1)gm(2) 
Or Ww) T Gn(Z) 


From Eqs. (10) and (19) the heat capacity above T, is then given by 


CAT) = Seams [0 + Il + Danla) - (+ aeoe) 


= [ø+ yr 2S - q ay? BHD) ya, (16) 


At T = T+, the chemical potential u = 0, so the above expression reduces to 


(rt) = [o+ patat - (y+ pe) ye (17) 


By comparing with Eq. (LJ), the discontinuity of heat capacity at the critical point is 


AC = O(T-) - C(T*) = re Ae yk (18) 


From the Eqs. (1), (17) and (18), we see that the heat capacity depends only on the particle 
number and the exponents p, l and q of the external potential. We can further examine the 
continuity of C (T) and its derivatives. As we know, g,(1) is divergent for 7 < 1. So, the heat 
capacity becomes discontinuous at the critical temperature T, when 7 > 1. From Eqs. (L1) 
and (16), the derivative 0C(T)/OT is discontinuous at T, for either 7 < 1 or 7 > 1. This 
is in agreement with Ref. [6]. From the analytical expressions we also see that the heat 
capacity at the critical point grows with increasing 7. 

For a rigid wall box, 7 = 1/2. It is straightforward to check that the above expressions 
indeed reduce to the known results in textbooks of statistical mechanics (e.g. Ref. [). As 


a special case let us consider the harmonic potential 


ie (19) 


We find the heat capacity C(T*), C(T7) and AC to be 


Cc 


O(T>) = oi (N)k, (20) 
a 1294(1) _ 993(1) 
eel | gs(1) - 20) Pi (21) 
— 993(1) 
Ac = A (Wyk. (22) 


This agrees with Refs. [2p]. 
Along similar lines we may analyze the thermodynamic properties of a trapped interact- 
ing Bose gas. In the mean field approximation, the single particle effective potential can be 


written as: 
Vestr) = V (r) + 2a7p(r), (23) 


where a is the s-wave scattering length, which characterizes the strength of interaction 
between two atoms. We assume |a| < A for the weak interaction considered here. By 
similarity with Eq. ([]), for the trapped interacting Bose gas the local grand potential can 


be written as 


pr) = —kTA* gsj2(€), (24) 


where 
€ = exp[—G(V(r) + 2ad*p(r) — p)]. (25) 


Expanding the function gs /2(€) in Eq. (24) according to the small parameter a/, and keeping 


terms only up to the lowest order in a/A, we have 


Ühele) ~ —kTA~ [g5/2(¢) — 2aX*p(r)g3,/2(6)] 
~ KTA | 95/2(6) — 2a" 95 0(6)] , (26) 
which coincides with the first order result of the virial expansion of the grand partition 


function [[3]. This is not surprising since the mean field approximation is accurate to the 


first order. From the grand potential we can obtain the local density 
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= 78 [g3/2(¢) — 4a g1/2(6) 93/2(0)] 


~ A~"gs/2(exp[—G(V(r) + 4a*p(r) — 1)]), (27) 


which, compared with the density of an ideal Bose gas, means the effective fugacity is 
exp|—6(V(r) + 4a\2p(r) — u)]. However, from a naive consideration based on the effective 
interaction one would think that it should be exp[—G(V(r) + 2aA?p(r) — u)| instead. The 
factor two in the atomic density expression naturally appears by taking the derivative of the 
grand potential with respect to u. 

As the counterpart to Eq. (B), we can find the total grand potential for the interacting 


gas at T > Te to be 


W Z 
Qpe = ~ igri [Gn42(2) — 2ad~* F372, 3/2, n(2)| , (28) 


where 
oo qiti 
= r (29) 
i,j=1 a J v(i F j)" i 
When a = 0, Eq. (28) recovers expression (H) for the non-interacting Bose gas, as it should. 


The average number of atoms is now given by 


= XG [Gms (2) = 2a~\F5/2,8/2n-1(2)| ; (30) 
where we have used 
OF: pala 1 
el) = 5 ena (31) 


As we know, g,+1(z) is a bounded, positive, monotonically increasing function of z, and 
F5y,(x) is positive. From Eq. BO) we see that the trapped interacting Bose gas has a higher 
chemical potential than an ideal Bose gas when the atomic number and temperature are the 


same in both systems. 


From Eq. (BO) for (N7) we can calculate the derivative 0u/0T at a fixed N‘, and the 


result is 
(4) _ E _ (nt 1)gmi(2) — 2aA7* (n + 8/2) F5/2,3/2,n-1(2) k 
oT (NI) T gn(z) — 2A} F3/2,3/2,n-2(2) 
H [ (n + 1)gn+1(2) 2a | (n + 1) Gn41(2) 3/2, 3/2, n-2(2) 
T In(2) À 9, (2) 
_ (n +3/2)F 32, 3/2, = | 3 (32) 
grz) 


Corresponding to Eq. (H), the total entropy of the trapped interacting Bose gas can then be 


derived as 


Wk 


= BA (7 + 2)gn+2(2) — BH gn+1(2)] 
2aWk 
igre lin + 5/2)F372, 3/2, n(2) — Bu F3/2, 3/2, 4-32)| (33) 


Calculating the derivative of the entropy ($/), we obtain the heat capacity above the BEC 


critical temperature to be 


i B os! 
CH(T)=T (= J 
2(2) 


E Ins 29n+1(2) 2a | (n + 3/2)(n + 5/2) F372, 3/2, n(2) 
= o PI a) =T at) Mesi A 
( 


+(n +1) jp _ eee] (N)k. (34) 


The last two formulas show the correction due to the mutual interaction of atoms to ther- 
modynamic quantities such as the entropy and heat capacity. 

The above derivations for the trapped interacting Bose gas is adequate only for the gas 
phase (without BEC), and furnish no information whatsoever about the condensed phase, 
or the nature of the condensation [[3]. Any attempts to use the same equations in the 
regime of BEC will lead to some mistakes. According to Ref. the mean field effective 
potential should be amended to include the contributions of atoms to Bose condensation. 


If we choose an appropriate effective mean field potential similar to Eq. (28) of Ref. [3], 
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we can, in principle, extend the above discussion to calculate thermodynamic quantities 
below the critical temperature. However, the derivation will be so involved that no simple 


analytical formulas can be expected to exist any more. 
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